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1. Introduction 



Phenomenology at the LHC often involves high multiplicity final states. For example, 
backgrounds to Higgs searches involve processes such as PP ^ W^W~ + 2 jets and 
PP ^ ti + bb. Both these examples involve 2 — 4 scatterings. At leading order (LO) 
such high multiplicity final state amplitudes can be evaluated using either numerical 
recursive techniques |jl|, || or other numerical and/or algebraic techniques 0, ^, H, 

However, O [as), next-to-leading order (NLO) corrections to the scattering am- 
plitudes are desirable. Not only do NLO corrections give a first reliable prediction 
of total rates, they also give a good error estimate on the shapes of distributions. At 
NLO the current state of the art for hadron colliders are 2 — > 3 processes. Thus NLO 
predictions for PP 3 jets |^ |10[ (based on virtual corrections of ref. p!3[| ) 
and PP —>■ V + 2 jets Jl^ (based on virtual corrections of ref. |jT5|, [16|, |1^) are 



known, and codes for PP ^ tt -\- jet |18, 19| and PP H + 2 jets via gluon 



fusion 1^ are under construction. Other processes such as PP V1V2 + jet and 
PP — > V1V2V3 are now feasible. 

By contrast the consideration of 2 — s> 4 processes is still in its infancy. In elec- 
troweak physics the full one-loop electroweak corrections to e~^e~ ^ 4 fermions were 
calculated in Ref. |2^ . However the calculation of NLO 2^4 QCD scattering 
cross sections is currently unexplored. Such a calculation involves both the evalu- 
ation of the one-loop six-point virtual corrections and the inclusion of the 2^5 
scattering bremsstrahlung contributions through Monte Carlo integration. 

In this paper we consider the virtual corrections to six-gluon scattering which 
is relevant for a calculation of PP 4 jets. By considering the one-loop correc- 
tions to gg gggg we select the most complicated QCD six-point processes. If the 
amplitude is calculated in terms of Feynman diagrams, the number of diagrams is 
very large and the gauge cancellations between these diagrams is the most severe. 
These cancellations could be a concern in a semi-numerical procedure; the six-gluon 
amplitude therefore provides a stringent test of the method. In this paper we con- 
sider neither the bremsstrahlung contributions, nor the one-loop processes involving 
external quarks, which are needed to obtain results for a physical cross section. 

The technique for the analytic calculation of the one-loop corrections to multi- 
gluon amplitudes which is relevant for this paper is the decomposition of the cal- 
culation into simpler pieces with internal loops of A/" = 4 and A/" = 1 multiplets of 
super-symmetric Yang-Mills particles and a residue involving only scalar particles 
in the loops [jTl], |2^. After recent advances p5| , |26| , p^ , p8| , all supersymmetric 



contributions have been computed analytically, however not all of the scalar contribu- 
tions for six-gluon amplitudes (or higher) are known yet. We present here numerical 
results for six-gluon contributions. For supersymmetric pieces we provide completely 
independent cross-checks of analytical results. 
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Although all one-loop 2^2 and almost all of the currently known 2 — > 3 
amplitudes were calculated using analytic techniques, we believe that semi-numerical 
or hybrid numerical/analytic techniques offer promise for more rapid progress. This 
technique was demonstrated recently for the case of the one-loop H -|- 4 partons 
amplitude [pO| . 

Many methods have been proposed to calculate NLO amplitudes, both semi- 
numerical [^, ^ |HT], ^ 0, ^ or numerical |D|, |T1. Of these 



methods only a few have actually been used to evaluate one-loop amplitudes. Only 
by using the methods in explicit calculations one can be sure that all numerical issues 
have been addressed properly. 

In section II we discuss the colour algebra involved with the evaluation of a six- 
gluon amphtude. The numerical techniques used in this paper are discussed in section 
III, while in section IV the comparison is made with numerous super-symmetric and 
the few scalar results, which exist in the literature. Finally, our conclusions in section 
V summarize the paper. 

2. Six-gluon amplitude at one- loop 

At tree-level, amplitudes with n external gluons can be decomposed into colour- 
ordered sub-amplitudes, multiplied by a trace of n colour matrices, T". The traceless, 
hermitian, Nc x Nc matrices, T", are the generators of the SU (Nc) algebra. Following 
the usual conventions for this branch of the QCD literature, they are normalized so 
that Tr(T''T^) = 5°^. Summing over all non-cyclic permutations the full amplitude 
A^j^^"^ is reconstructed from the sub-amplitudes A^^'^'^{a) 



ATi{p^^^,<^^})=9''~' E Tr (T^D . . . TM™, ) (^'(1^ • • • , p'^)' ) • (2.1) 

n 

The momentum, helicity (±), and colour index of the z-th external gluon are denoted 
by Pi, Aj, and a, respectively, g is the coupling constant, and Sn/Zn is the set of 
{n — 1)\ non-cyclic permutations of {1, ... , n}. 

The expansion in colour sub-amplitudes is slightly more complicated at one- 
loop level. Let us consider the case of massless internal particles of spin J = 0, 1/2, 1 
corresponding to a complex scalar, a Weyl fermion or a gluon. If all internal particles 
belong to the adjoint representation of SU(Ai'c), the colour decomposition for one-loop 



n-gluon amplitudes is given by P3 



Ln/2J+1 

A^^{{v^,Ka,})=g- E Gr„»4;I(^), (2.2) 

C=l (tGSu/ Sn;c 

where [x\ denotes the largest integer less than or equal to x and Sn-c is the subset of 
Sn which leaves the double trace structure in Gr„.c(l) invariant. The leading-colour 
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structure is simply given by, 



Gr„;i(l) = Nc Tr(T'^i ■ • ■T'^") . (2.3) 
The subleading-colour structures are given by products of colour traces 

GinM) = Tr(T"i ■ ■ ■ T"=-') Tr(T"= ■ ■ ■ T'^") . (2.4) 
The sub leading sub-amplitudes An-c>i are determined by the leading ones A^^^'j^ 



through the merging relation |2^, ^ 



<Ui(l,2,...,c-l;c,c+l,...,n) = i-iy-' AL;'i(ai,...,a„), (2.5) 

aeOP{a}{f3} 

where a, G {a} = {c — 1, c — 2, . . . , 2, 1}, Pi G {/?} = {c, c + 1, . . . , ri — 1, n}, 
and 0P{a}{/3} is the set of ordered permutations of {1,2,..., n} but with the last 
element n fixed. The ordered permutations are defined as a set of all mergings of at 
with respect to the Pi, such that the cyclic ordering of the aj within the set {a} and 
of the Pi within the set {P} is unchanged. In practice, since n is fixed, no further 
cycling of the set {P} is required. Thus a complete description can be given in terms 
of the leading colour sub-amplitudes An-i alone. 

The contribution of a single flavour of Dirac fermion in the fundamental repre- 
sentation, (relevant for quarks in QCD) is 

A°'^^^({p.,A.,a,})=^" J2 MT''^^'^---T''^^-^)A%'\pll^\...,pti::;). (2.6) 

cr£S„/Z„ 



Simple colour arguments allow one to demonstrate that this colour sub-amplitude 
is the same as the leading colour sub-amplitude for a single Weyl fermion in the 
adjoint representation defined in Eq. (2.2). 

Since the sub leading colour amphtudes are not independent, we shall henceforth 
drop them from our discussion. To simplify the notation we shall also drop the 
subscripts n and c. The amplitude denoted by A will thus refer to leading colour 
amplitude with six external gluons. 

3. Method of calculation 

The method we use is purposely kept as simple as possible. Especially in numerical 
methods this is desirable for both keeping track of numerical accuracy and code 
transparency. 

To generate all the required Feynman diagrams we use Qgraf |4^. The Qgraf 



output is easily manipulated using Form to write the amplitude in the form 



N 



A{1, 2, 3, 4,5,6) = J2Y1 ^/^--/^M iPu ei; • • • ; Pe, e,)!^^'-^"^' {p,, . . . ,p,) , (3.1) 



N=2 M=0 



3 



where the kinematic tensor K depends on the purely four-dimensional external vec- 
tors and contains all the particle and process information. The A^-point tensor inte- 
grals of rank M are defined in D dimensions as 



,P6) 



in^/'^ did2- ■ - d 



N 



di = {1 + Qi 



(3.2) 



and can be evaluated semi-numerically. 

For < 4 we use the method of ^ which we already developed, 

tested and used in the calculation of H -|- 4 partons at one-loop In general, 

the basis integrals will contain divergences in e = {4 — D)/2 from soft, coUinear and 
ultraviolet divergences and the answer returned by the semi-numerical procedure will 
be a Laurent series in inverse powers of e. 

For the five (six)-point tensor integrals the method we use relies on the complete- 
ness (over-completeness) of the basis of external momenta for a generic phase space 
point. We therefore use a technique for tensor reduction which generalizes the meth- 
ods of ref. [0, ^]. This technique is valid as long as the basis of external momenta 
is complete^. Assuming we have a complete basis of external momenta we can select 
a set of 4 momenta {pfci , Pfc2 ? Pfca , Pfc4 } which form the basis of the four- dimensional 
space. We can then decompose the loop momentum 



4 4 



(3.3) 



i=l 



i=l 



where the f fc. are defined as linear combinations of the basis vectors 



where G is the Gram matrix and 

V = -^J2^rk^-rk^_i)v',^^, rk = ql 

i=l 



(3.4) 



(3.5) 



With this relation it is now easy to reduce an A^-point function of rank M to a lower 
rank A^-point function and a set of lower rank (A^ — l)-point functions 



T^l^■■■^lM 



1 ^ 

9 E (^"^^ 



') 



(3.6) 



^For exceptional momentum configurations (such as threshold regions or planar event configu- 
rations) this is not the case. Exceptional configurations can be treated using a generalization of the 
expanded relations proposed in refs. Esl Eol. This is beyond the scope of this paper. 
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where Inj is a (A^— l)-point integral originating from In with propagator dj removed. 
More exphcitly, choosing without loss of generality the base set {^1,^2,^3,^4}, we 
get 



^j^'"^''*(Pl,P2,P3,P4,P5, ■■■,Pn) = In ^ \Pl,P2,P3,P4,P5, ■ ■ ■ , Pn)V'"" {pi, P2, Ps, p^) 

1 
2 



+ l (^TV-T' '{Pl+P2,P3,P4,P5, ...,Pn)- In-T^" '{P2,P3,P4,P5, ■ ■ ■ ,Pn)) 



X<''(pi,P2,P3,P4) 
1 

2 



+ J (^TV-T'' ' (Pl' P2 + P3, P4, P5, • • • , Ptv) - In-T"' ' {Pl + P2,P3, Pi, Ps, • • • , Pn)) 



1 

2 



+ ^ {1""' ' (Pl, P2, P3 + P4, P5, • • • , Ptv) - /j^™' ' {Pl,P2 + P3, P4, P5, • • • , Pn)) 



^vTi.Pl,P2,P3,Pi) 

+ o (^TV-T'"' {PUP2,P3, P4+P5,- ■ - ^Pn) - 1'^'^'^^'"'^ {PUP2,P3 + Pi, PS, • • • , Pn)) 



2 

XWr'(pi,P2,P3,P4) • 



(3.7) 



For example, applying this relation repeatedly to the tensor six-point integrals 
we will be left with the scalar six-point integral and five-point tensor integrals. The 
five-point tensor integrals can be reduced using the same technique. Subsequently 
we can use the method of ^ ^ to further numerically reduce all remaining 



integrals to the basis of scalar 2-, 3- and 4-point integrals. This procedure turns out 
to be efficient and straightforward to implement numerically. 

4. Comparison with the literature 

Since we have directly calculated the loop amplitudes with internal gluons and 
fermions we can easily obtain the result for QCD with an arbitrary number rif of 
flavours of quarks, 

^QCD ^ ^[1] ^ ^^[1/2] _ 

However since the analytic calculations in the literature are presented in terms of 
supersymmetric theories we need to re-organize our results to compare with other 
authors. 

4.1 Supersymmetry 

Since we have calculated the amplitudes with massless spin 1, spin 1/2 and spin 
particles in the internal loop we can combine our results as follows 

A^=4 = A[i]+4yl[i/2] + 3^[o]^ (4.2) 
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A^^'^, so constructed, describes an amplitude where the full supersymmetric Af = A 
multiplet runs in the loop, and A'^^^ denotes the contribution from an A/" = 1 
super-multiplet running in the loop. 

In analytic calculations the intention is to proceed in the opposite direction. 
Amplitudes with multiplets of supersymmetric Yang-Mills in internal loops have 
much improved ultra-violet behavior and are four-dimensional cut-constructible. For 
this reason, all of these supersymmetric amplitudes have been calculated and most 
have been presented in a form suitable for numerical evaluation. As far as six-gluon 
amplitudes with scalars in the loop, A'^^\ are concerned three of the needed eight 
independent helicity amplitudes have been published so far. Only in the helicity 
combinations where all contributions are known can one reconstruct the ingredients 
needed for QCD amplitudes 

^[1] = A^=' - 4A^=' + AM , (4.4) 

4.2 Numerical results 

As a preparatory exercise we performed a check of the four- and five-point gluon 
one-loop amplitudes. We found agreement with the literature [|^, |5^, [Tl[ . 

We now turn to the amplitude for six-gluons which is the main result of this 
paper. Our numerical program allows the evaluation of the one-loop amplitude at 
an arbitrary phase space point and for arbitrary helicities. For a general phase space 
point it is useful to re-scale all momenta so that the momenta of the gluons, (and the 
elements of the Gram matrix), are of 0(1) before performing the tensor reduction. 
Without loss of generality we can assume that this has been done. 

To present our numerical results we choose a particular phase space point with 
the six momenta pi chosen as follows, {E,px,Py,Pz), 

Pi = ^(— 1, -|- sin^, + cos ^ sin -|- cos6' cos</)), 

P2 = ^(— 1, — sin^^, — cos6'sin0, — cos6'cos0), 

P3 = |(1,1,0,0), 

P4 = ^(l,cos/5,sin/5,0), 

P5 = —(1, COS a COS /5, COS a sin /5, sin a), 
6 

Pe = -Pi - P2 - P3 - P4 - P5 , (4.6) 

where 9 = 7r/4,0 = 7i/6,a = 7r/3,cos/3 = —7/19. Note that the energies of pi and 
P2 are negative and pf = 0. In order to have energies of 0(1) we make the choice for 
the scale fi = n = 6 [GeV]. As usual /i also denotes the scale which is used to carry 
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the dimensionality of the D-dimensional integrals. The results presented contain no 
ultraviolet renormalization. 

Analytic results require the specification of eight helicity combinations: all other 
amplitudes can be obtained by the parity operation or cyclic permutations. We 
choose these eight combinations to be the two finite amplitudes (+ + + + ++, — + + + 

++), the maximal hehcity violating amplitudes ( — I — hH 

and the next-to- maximal hehcity violating amplitudes ( 1 — I- 

— I h). These eight amplitudes would not be sufficient for a numerical evaluation, 

but the numerical approach allows the evaluation of any helicity configuration at 
will. 

In Table |l]we give results for a particular colour sub-amplitude A'^^^{1^ 2, 3, 4, 5, 6) 
for the above eight choices of the helicity. An overall factor of icr has been removed 
from all the results in the Tables |1|, and ^ 

_ (47rrr(l + e)r^(l-6) 
^ 167r2 r(l-2e) ' ^ ' 

The results for the A/" = 4 amplitudes depend on the number of helicities of gluons 
circulating in internal loops. For a recent description of regularization schemes see, 
for example, ref. |Q. Our results are presented in the 't Hooft-Veltman scheme. The 
translation to the four-dimensional helicity scheme is immediate 

^FDH = ^t^HV + "^^tree • (4.8) 

Note that analytic results from the literature are quoted in the four- dimensional 
helicity scheme, which respects supersymmetry. These results have been translated 
to the 't Hooft-Veltman scheme using Eq. ( [4. 81) before insertion in our tables. 

In Table |^ we give results for the colour sub-amplitudes ^•'^^^(l, 2, 3, 4, 5, 6) for 
the same eight helicity choices and where possible compare with analytical results. ^ 
Note that because of the relation 

-A Isingular ^ ^ ; (4-9) 

the column giving the single pole can as well be considered as a listing of the results 
for the colour-ordered sub-amplitudes at tree graph level (stripped only of the overall 
factor of z). 

We note that for two of the helicity amplitudes 1 h+ and — I 1 h 

we were unable to evaluate the analytic results numerically. This was due to the fact 
that calculating the residue of certain poles as required by the formula in ref. [p^ ], 
resulted in zero value denominators of sub-expressions^. 

^In Eq. (5.16) of ref. |Q for the degenerate case ni=j-l=2 one has Cm = {j + ■ ■ ■ ,n — 1}, as 
can be seen from Fig. 8 of this same paper. This point has also been made in ref. ||55|] . 

■^We thank the authors of ref. ||2^ for confirming that there are problems with the numerical 
evaluation of the formula for these amplitudes in their paper. 
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Hclicity 




1/^ 


1 


[Rcf]/(Eq.#) 


+ + + + ++ 
++++++ 




(-1.034 + i 2.790)10"* 




(-9.615 + i 3.708)10"* 




-(0.826 + i 2.514)10"^ 


[SN-A] 


h + + ++ 

-+++++ 




(1.568 + j 2.438)10"** 




(-0.511 + i 1.129)10"^ 




-(3.073 + i 0.1223)10"^ 


[SN-A] 


h + ++ 

h + ++ 


-161.917 + i 54.826 
(-0.933 + j 1.513)10"* 


—489.024 — i 212.415 
-(7.655 + i 0.440)10"* 


—435.281 — i 1162.971 
-(-0.221 + j 1.834)10"^ 


|23)/(4.19) 

[SN-A] 


- H h ++ 

- H h ++ 


—33.024 + i 44.423 
{-7.542 + i 0.939)10"* 


-169.358 + i 33.499 
-(1.157 + i 0.363)10"* 


—330.119 — i 229.549 
-(3.474 + j 2.856)10"* 


|23)/(4.19) 

[SN-A] 


h H h + 

- + H h+ 


—0 5720 — ?■ 3 939 
(-2.279 + j 1.803)10"* 


fi qoQ — )■ 1 n 302 
-(1.176 + i 0.399)10"^ 


28 469 — 7 5 058 
(0.054- j 3.307)10"^^ 


E3I l(A 19") 
[SN-A] 




h ++ 


-6.478 - i 10.407 
(2.686 - i 1.668)10"* 


6.825 - i 37.620 
(1.232 + i 0.554)10"^ 


75.857 - i 47.081 
(0.020 + i 3.334)10"^^ 


|24)/{6.19) 

[SN-A] 


h - ++ 

h - ++ 


14.074 - i 22.908 
-(1.619 + i 0.943)10"* 


80.503 - i 23.464 
-(1.030 + j 8.234)10"* 


169.047 + j 93.601 
(1.560 -i 0.801)10"* 


j24|/{6.24) 

[SN-A] 


- H h - + 

- H h - + 


13.454 + j 13.177 
(1.045 - i 0.113)10"^ 


3.495 + i 58.632 
(-0.772 + i 1.652))10"* 


-88.32 + i 103.340 
(-7.795 + i 7.881))10"* 


|24|/{6.26) 

[SN-A] 



Table 1: AA=4 color ordered sub- amplitudes evaluated at the specific point, Eq. (^). 
The results are given in the 'tHooft-Veltman regularization scheme. [SN-A] means the 
difference between the semi-numerical result and the analytical one. 



Hclicity 




lA 


1 


[Ref]/(Eq.#) 


++++++ 
++++++ 




(-3.470 + i 9.320)10"^ 




(-3.226 + j 1.253)10"* 




-(3.899 + i 8.969)10"* 


[SN-A] 


-+++++ 
-+++++ 




(5.228 + j 8.127)10"^ 




(-1.678 + 1 3.775)10"* 




-(1.013 + i 0.2066)10"'^ 


[SN-A] 


h + ++ 

h + ++ 




(-3.297 + i 5.194)10"^ 


26.986 -j 9.1376 
-(-2.104 + i 0.344)10"* 


101.825 -i 52.222 
(0.949 - i 4.895)10"* 


0/(5.9) 
[SN-A] 


-+-+++ 
-+-+++ 




(-1.847 + 1 0.8566)10"!" 


5.504 - i 7 Am 
-(6.141 + j 4.633)10"!° 


21.811 - i 29.051 
(3.095 + j 2.138)10"^^ 


13/(5.12) 
[SN-A] 


-++-++ 
-++-++ 




(-7.599 + j 6.018)10"^ 


0.09533 + i 0.6565 
-(3.929 + i 1.304)10"* 


-2.183 + i 3.260 
(0.008 - i 1.100)10"^^ 


0/(5.12) 
[SN-A] 


+ ++ 

+ ++ 




(8.965 -i 5.555)10"^ 


1.080 + i 1.735 
(4.107 + i 1.858)10"* 


0.722 + i 5.285 
(0.002 + j 1.114)10"^^ 


[§/(9) 
[SN-A] 


\ h+ 

\ h+ 




(-5.351 - i 2.825)10"^ 


-2.346 + i 3.819 
-2.346 + i 3.819 


-2.238 + i 17.687 


11/(5.4,2.3) 

[SN] 


- H 1 h 

- + - H h 




(1.124 - i 0.2060)10"!° 


-2.242 - i 2.196 
-2.242 - i 2.196 


-1.721 - i 7.433 


11/(5.13,2.3) 

[SN] 



Table 2: M=\ color ordered sub-amplitudes evaluated at the specific point, Eq. ( [4. 61) . [SN] 
means that the result is obtained using our semi-numerical code, while [SN-A] denotes the 
difference between the semi-numerical result and the analytical one. 



Lastly in Table |^ we give results for the colour sub- amplitudes A'^l (1, 2, 3, 4, 5, 6) 
for scalar gluons, for the same eight helicity choices.^ For all amplitudes for which no 
analytic result exists, we checked the gauge invariance of the ampUtudes by changing 
the gluon polarization. The gauge invariance was obeyed with a numerical accuracy 
of O (10^*^). To evaluate a single colour-ordered sub-amplitude for a complex scalar 
took 9 seconds on a 2.8GHz Pentium processor. To evaluate the complete set of 

^In ref. p7[ ] [vl-v3] the definition of Ff has an overall sign missing, a typographical error not 
present in the original calculation of the M — 1 term in ref. [e3 . 
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Hclicity 


1 / 2 




1 


[Rcf 


/(Eq.#) 


III I J 1 

i — \ — 1 — 1 — \ — r 

+ + + + ++ 


n 
U 

(3.672 + i 9.749)10"^ 


u 

(— 3.404 + i 1.238)10-8 


— (3.016 + i 9.169)10-8 


§ 

JS 


/(4.3) 
N-A] 


h + + ++ 

-+++++ 




(5.921 + i 8.411)10"^ 




(-1.606 + i 4.051)10-8 


—3.194 + i 0.6503 
-(1.086 + i 0.038)10-^ 


|26|/{4.10) 

[SN-A] 


h + ++ 

--++++ 




(1.280 + i 0.002)10-8 


8.995 - i 3.046 
(2.768 + i 4.232)10-8 


43.089 - i 20.288 
(-1.004 + j 0.955)10-^ 


1^/(4.27,4.28) 

[SN-A] 


- H h ++ 


(1.045 - i 0.580)10-8 


1.835 - i 2.468 


9.752 - i 11.791 


[SN] 


-++-++ 


(-7.791 + 1 6.717)10-^ 


3.178 ■ 10-^ +i 0.2188 


-1.447 + i 0.1955 


[SN] 


h ++ 


(8.934 - i 5.359)10-s 


0.3599 + i 0.5782 


0.5617 + j 5.8166 


[SN] 


\ h+ 


(0.1016 + j 1.276)10-8 


-0.7819 + i 1.273 


-0.6249 + i 6.552 


[SN] 


- H h - + 


(1.065 - i 0.5417)10-8 


-0.7475 - i 0.7321 


-1.298 - i 3.255 


[SN] 



Table 3: One loop six gluon colour ordered sub-amplitudes with a scalar loop evaluated 
the specific point Eq. ([4.6D . [SN] means that the result is obtained using our semi- numerical 
code, while [SN-A] denotes the difference between the semi-numerical result and the ana- 
lytical one. 

64 possible helicities will be less than 64 times longer, because the scalar integrals 
are stored during the calculation of the first amplitude are applicable to all other 
configurations with the same external momenta. 

5. Conclusions 

In this paper we have presented numerical results which demonstrate that the com- 
plete one-loop amplitude for six-gluon scattering is now known numerically. By 
forming multiplets of SUSY Yang Mills in the internal loops, we were able compare 
with most of the known analytic results. In addition, we have presented numerical 
results for amplitudes which are currently completely unknown. Note that the ana- 
lytic and semi-numerical results are complementary. The hardest piece to calculate 
analytically is the scalar contribution A^, which is the easiest for the semi-numerical 
approach. Thus it is possible that a numerical code involving both semi-numerical 
and analytic results will be the most efficient and expedient. Our results demonstrate 
the power of the semi-numerical method, which can supplant the analytic method 
where it is too arduous and provide a completely independent check where analytic 
results already exist. 

After inclusion of the one-loop corrections to the other parton subprocesses in- 
volving quarks it would be possible to proceed to a NLO evaluation of the rate for 
four jet production. We intend to use these methods to calculate NLO corrections to 
other processes which we consider to be of more pressing phenomenological interest. 
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